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Abstract 

The thermally activated flips of the local spontaneous polarization in relaxors 

were simulated to investigate the effects of the applied-ac-field amplitude on 

the dielectric susceptibility. It was observed that the susceptibility increases 

with increasing the amplitude at low temperatures. At high temperatures, 

the susceptibility experiences a plateau and then drops. The maximum in the 

temperature dependence of susceptibility shifts to lower temperatures when 

the amplitude increases. A similarity was found between the effects of the 

amplitude and frequency on the susceptibility. 
PACS: 77.22.Ch, 77.80.-e, 77.84.-s, 77.84.Lf 
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Relaxor ferroelectrics (relaxors) have been studied for nearly 40 years since 
Pb(Mg 1/ / 3 Nb 2 / 3 )03 (PMN) was first synthesized by Smolenski et a/..0The dielectric response 
of relaxors is characterized by the diffuse phase transition (DPT) and a strong frequency 
dispersion.il Various models, such as the compositional heterogeneity model,0 the superpara- 
electric model,! and the glasslike model J! et al, were proposed to rationalize the complicated 
behaviors of relaxors. It is widely accepted nowadays that the presence of polar microregions 
m nanosca is responsible for the relaxor behaviors. 

The effects of the applied ac field on relaxorsEH3 cause great interest since they provide 
some clue of the relaxation mechanism. Glazounov et al. observed that the dielectric per- 
mittivity of PMN increases with increasing amplitude of the applied ac field.B A similarity 
was also found between the effects of the amplitude and frequency on the permittivity. In 
addition, the ac-drive-enhanced relaxor characteristics and domain breakdown were observed 
in (PbLa)(ZrTi) (PLZT).0 There are two possible mechanisms, i. e., domain- wall motion 
model and superparaelectric model, to explicate the nonlinearity of dielectric permitivity of 
PMN relaxors. Glazounov et al. sugg estedi that it is related to domain-type process rather 
than thermally activated flips of the local spontaneous polarization (i.e. superparaelectric 
model). However, they did not consider the interaction of polar microregions when investi- 
gating the superparaelectric model, which is just one of the key points related to response 
of the external fieldillHl!. In this study, we conduct a Monte Carlo simulation to investigate 
the influence of measuring field on the dielectric susceptibility of relaxors. 

We investigate the thermally activated flipping process of the local spontaneous polar- 
ization in relaxors. Following the work of Gui et a/.,0 the polar microregions are regarded 
as point dipoles. Then relaxors are modeled to be a system consisting of Ising-like dipoles 
with randomly distributed interactions:^ 

H = Jij OiOj - E ext J[J2 ^' C " S — (1) 

i# i V 

where a,, aj = ±1 are dipole spins. When the projection of the ith dipole moment /I, on 

the direction of the external field E ext is positive, cij takes value +1, otherwise, <Ji takes 



value -1. Oi is the angle between /Zj and E ext , and fx is the maximal magnitude of the dipole 
moments. is the effective interaction energy between the nearest neighbor dipoles, which 
has a Gaussian distribution with a width AJ. reflects the correlation between polar 
microregions, which is essential to the glassy behaviors. In general, the external field 
contains a measuring ac field and a bias dc field. In this paper, only the ac field is involved, 
i.e., 

E ext = E exp (* 27r ^-) > ( 2 ) 

where t is the real time. E and are the amplitude and the period of the ac field, 
respectively. 

The Monte Carlo simulation is performed on a 16 x 16 x 16 simple cubic lattice with 
periodic boundary conditions. The details of simulation process can be found in Ref. 15. 
The dielectric susceptibility is defined as 

t{)-\-t i, s . , 

/tts / P(0exp(i27r£)dt. 
X = C{ * , (3) 

where C is a proportional factor which is chosen to be 1 in this contribution, and (• ■ ■) 
denotes the configurational averaging. p(t) is the normalized polarization: 

i ^ 

During the simulation process, p(t) is recorded and x is calculated according to Eq. (3). 
We choose to = 200MCS/dipole to eliminate the influence of the initial state and t b s = 
3000MCS/dipole to be the observation time. The simulation is performed in many runs 
with different initial conditions so that the configurational averaging can be done. Longer 
observation time was also adopted in test, but no obvious influence on results was observed. 

In order to verify the validity of the method, the dielectric susceptibility under a weak 
field is firstly calculated. The result is shown in Fig. 1. It can be seen that the susceptibility x 
reaches its maximum at a certain temperature (T m ) and changes gradually around T m , which 
is known as the diffuse phase transition (DPT) in relaxors. A strong frequency dispersion can 



be also observed: x decreases with increasing field frequency at low temperatures, and T m 
moves to higher temperatures. All these characteristics are consistent with the experiments! 
and the previous theoretical results^. 

Now, let us investigate the effects of the field amplitude on the dielectric susceptibility. 
The susceptibility curves under different ac-field amplitudes Eq are depicted in Fig. 2 when 
the measuring frequency is kept as ti = lOMCS/dipole. (We express the frequency by 
here and hereafter.) From Fig. 2 one can list the most essential features of the nonlinear 
effect: (1) the dielectric susceptibility increases with increasing Eq at temperatures T < T m 
where the frequency dispersion is observed; (2) increasing Eq will make the maximum in the 
temperature dependence of x shifts to lower temperatures, which has the similar effect of 
decreasing frequency (see also Fig. 1). The change of the imaginary part, x" i shows similar 
features in the simulation. These features agree with the experiments in pmnII very well. 
The concepts of "slow dipole" and "fast dipole" can help to understand the increasing of 
the susceptibility. Slow dipoles are those dipoles which flip too slow to keep up with the 
changing of the ac field and give no or little contribution to the dielectric susceptibility. At 
low temperatures, there are large amounts of slow dipolesJil When Eq increases, the driving 
force on slow dipoles is enhanced. Slow dipoles are forced to flip faster and they give more 
contribution to the dielectric susceptibility x- F° r fast dipoles, the contribution changes 
slightly at low drives (see below). As a result, the susceptibility \ increases with increasing 

Eq. 

It can be seen in Fig. 2 that the dielectric susceptibility slightly decreases with in- 
creasing the external-field amplitude Eq at high temperatures. The tendency is weakened 
at higher frequencies while becomes more evident at lower frequencies. Fig. 3 shows the 
cases for a lower frequency t L = 50MCS/dipole. It shows that the dielectric susceptibility 
decreases at high temperatures and increases at low temperatures when Eq increases. These 
results are similar to the experimental cases in PLZlfflto some extent. However, the com- 
puted maximum in x{T) decreases with increasing Eq, which is opposite to the experimental 
observations.0 It reflects the defect of the model or/and the method we used. 



To get further knowledge of the continuous effects of the ac-field amplitude, we plot in 
Fig. 4 the curves of \ as functions of the amplitude E for different temperatures when the 
measuring frequency is fixed as = lOMCS/dipole. At low temperatures, the dielectric 
susceptibility increases first, and then drops with increasing Eq. This means that the applied 
field speeds up the flipping of dipoles at small E values so x increases first, while the system 
is nearly saturated at large Eq values which causes the drop of x- At high temperatures, the 
dielectric susceptibility experiences a plateau at the beginning and then decreases when the 
applied field increases. These results are consistent with the experiments in PMN when Eq 
varied in wide range of values.0 In Ref. 7 and Ref. 8, Eq is not large enough, so x increases 
at low temperatures and remains steady at high temperatures with increasing E . 

Fig. 5 demonstrates the field dependence of x & t different measuring frequencies and a 
fixed temperature T = 1.5AJ/&B. It shows that the maximum of the curve shifts to lower 
field amplitude when decreasing the measuring frequency. The shapes of curves are similar 
for different frequencies. 

Fig. 6 shows the temperature of the susceptibility maximum (T m ) as a function of the 
external- ac-field amplitude Eq. A nonlinear relation can be found between T m and Eq. It is 
conflict with the linear law observed in experiments^!!. Perhaps the field used in experiments 
is not large enough to reveal the high-order effects of the T m ~ Eq curve. Further experiments 
are needed to testify the theoretical predictions. 

The Eq. (3) could be generalized to include the Fourier component at different frequencies 
than that of E ext . Fig. 7 gives the curve of X2^IXi where X2u is the second-order component 
of the susceptibility. It can be seen that X2^/x is stronger at lower E ext and T. 

By means of the results above, we can see that the behaviors of the system described by 
the model Hamiltonian in Eq. (1) are consistent with many aspects of the experiments when 
the applied-ac-field amplitude varies. There are two points that should be mentioned here. 
First, the interactions between polar microregions play an important role in the dielectric 
response. If the interaction does not exist, the dielectric susceptibility will decrease with 
increasing field amplitude as what is pointed out by Glazounov et a/.Bfi Secondly, the model 
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in Eq. (1) is a rather simplified model. It cannot reflect the effects of external field on the 
crystal structure comple tely.0 It describes only the thermally activated flips process of the 
local polarization. Indeed, there may be more dielectric mechanism in relaxors. For example, 
It was presented that there may be two kinds of polarization processes in relaxors.0 Very 
recently, various types of contributions were found to dominate the dielectric response within 
different ac-drive amplitude rangesJll 

In conclusion, the simulation results suggest that the thermally activated flips of the local 
spontaneous polarization in relaxors plays an important role in producing the relaxation 
phenomena. 

This work was supported by the Chinese National Science Foundation. 
* Electronic address: zrliu@phys.tsinghua.edu.cn 
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FIGURES 

FIG. 1. Weak-field susceptibility as a function of temperature (in units of AJ/ks)- The 
field amplitude is fixed as £"0=0.1 AJ//Z. The curves 1-4 correspond to the field frequency 
t L = 100, 50, 20, lOMCS/dipole, respectively. 

FIG. 2. Dielectric susceptibility at various field amplitudes, E (1-0.5, 2-1.0, 3-1.5, 4-2.0AJ//2). 
The field frequency is kept as tL = lOMCS/dipole. Inserted graphics is the imaginary part of 
susceptibility. 

FIG. 3. Dielectric susceptibility at various field amplitudes, Eq (1-0.5, 2-1.0, 3-1.5, 4-2.0AJ//i). 
The field frequency is kept as ti = 50MCS/dipole. 

FIG. 4. Field amplitude (in units of AJ//2) dependence of dielectric susceptibility at a fixed 
field frequency ti = lOMCS/dipole. Curves 1-5 correspond to temperatures T=0.5,1.5,2.5,3.5 and 
4.5A J/ks, respectively. 

FIG. 5. Field amplitude dependence of susceptibility at various frequencies and a fixed tem- 
perature T = l.5AJ/ks- The amplitude and the frequency are measured in units of AJ/fi and 
MCS/dipole, respectively. 

FIG. 6. Temperature T m , corresponding to the position of the maximum in x(T), as a function 
of the field amplitude (in units of AJ/jl). The temperature is measured in units of AJ/ks- 

FIG. 7. Temperature dependence of X2w/x at various field amplitudes, Eq(1-0.5,2-1.5AJ/ jl) . 
The field frequency is kept as ti,=10MCS/dipole. 
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